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A quantum spin system can be modelled by an equivalent classical system, with an effective 
Hamiltonian obtained by integrating all non-zero frequency modes out of the path integral. 
The effective Hamiltonian H e g ({S;}) derived from the coherent-state integral is highly sin- 
gular: the quasiprobability density exp(— f3H c f{ ), a Wigner function, imposes quantisation 
through derivatives of delta functions. This quasiprobability is the distribution of the time- 
averaged lower symbol of the spin in the coherent- state integral. We relate the quantum Monte 
Carlo minus-sign problem to the non-positivity of this quasiprobability, both analytically and 
by Monte Carlo integration. 



1 Introduction 

It is possible, therefore, that a closer study of the relation of classical and 
quantum theazry might involve us in negative probabilities, and so it does. 
R P Feynmanlil 

One may distinguish two broad numerical approaches to quantum statistical me- 
chanics. One can evaluate a path integral by direct Monte Carlo methods. Here 
the notorious minus-sign problem (or phase problem) often hinders direct evalua- 
tion of the path integral: this refers to a rapid oscillation of the integrand in sign or 
phase, and the resulting intolerably slow convergence of integrals evaluated by ran- 
dom sampling. Alternatively, one can integrate out the quantum fluctuations, leaving 
an effective. Hamiltonian H e g(x) with c-number variables x to be used in a classical 
simulation^. A similar non-positivity arises here if the observables described by x 
are incompatible: the distribution exp(— {3H c g (x)) is then a Wigner function, which 
is not in general positive-definite. 

A free spin s serves here as a useful toy model to relate the path integral and 
the Wigner functions. The spin coherent-state path integral, a poor starting point for 
numerical simulations, presents the worst case of the sign problem: for continuous 
paths the integrand is a pure Berry phase factor. The spin Wigner function is highly 
singular, consisting_nf derivatives of delta functions supported on concentric spheres 
of quantised radiustm We demonstrate the common origin of these non-positivities 
both analytically and numerically. 
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2 Classical effective Hamiltonian 



We consider a system with density matrix p = exp(—0H), and assign c-number 
variables x £ M. N to operators x. We define the Wigner function W{x) and classical 
effective Hamiltonian H g (x) of these variables as 

W(x) = exp(-/3if eff (x)) = Tr (p5 N (x - x)) , (1) 

where Sn is the TV-dimensional delta function. To remove operator-ordering ambi- 
guity, we define 



>N 



f d A 

(x-x) =J j^yT exp(iA- (x-x)). (2) 



If the components of x do not commute, the delta operator need not be positive. 

Here we restrict consideration to the spin operator x = S in the spin-s represen- 
tation. The distribution is the spin Wigner function W S (S). Ear a single spin with 
vanishing Hamiltonian the trace of Eq. (||) is easily evaluatedBa to give derivatives 
of delta functions supported on concentric spheres of quantised radius: 



(3) 



where S = |S| and S 2 = s(s + 1). The spheres with m < do not contribute and 
for integer spin there is a positive delta function at the origin of weight 

-S'(S) S 3 (S) (4) 



2ttS(2s + 1) (2s + 1) 

To motivate this form, we need to verify that the correct marginal distribu- 
tions are obtained. Integrating W s over a plane yields the exnected distribution 
12m=-s ' ^ — m)/(2s + 1) for the normal component of spind. For two coupled 
spins 1 /2, with Hamiltonian — JSi • S2, a convolution integral of the two single-spin 
functions W1/2 gives (after some manipulations) 

W{b 1 + b 2 )- 3e/ 3.//4 + e -3/3J/4 • ^ 

This is the correct thermal average of triplet and singlet distributions, in line with 
the effective-Hamiltonian interpretation. Eq. (Q) can be further interpreted as a lo- 
cal hidden variable distribution, which must give the correct correlations for spin 
measiirements on two electrons. No positive definite distribution exists for a singlet 
statea Q, and the distribution must therefore be non-definite. 
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3 Coherent-state path integrals 



Spin coherent states |n) for spin s are labelled by a unit vector n, such that 
n • S|n) = s\n). To relate the Wigner functions just described to the coherent state 
path integral, we split the exponent of Eq. (||) into L time slices, 



W S (S) = (5 3 (S - S)> = J 0e' A S / (l - iX ■ S/L + 0(L" a )J ) . (6) 

We now insert a resolution of unity in the form ^§g=r J dn\n) (n| between each slice!, 
represent the spin by its lower symbol (s + l)rJ3'El 

S = ^^ /rfn( S + l)n|n)(n|, (7) 

and re-exponentiate. The lower symbol is used as the paths (here and in our Monte 
Carlo calculation) are only piecewise continuous. This gives 

W a (S) = Urn W S L (S), (8) 

L — >oc 

the limit defined in the sense of a distribution, with approximants 

-i(s+l)A-n 



. L 

2 



3 x ' 2s + 1 / 8tt 3 / 11 V 4?r 

1=1 v 

X (m | n 2 ) (n 2 \---\n L ) (n L \ n x ) (9) 



2s + 1 

where fi = Ya=i n i/L is the time-averaged coherent-state label. We therefore have 
a sequence of well-defined distributions (g) convergent on the Wigner function. For 
large L the discretised distribution W a (S) can be shown to be asymptotically a 
smeared Wigner function; for spin 1/2 we find 

</ 2 (S) « (^) 3/2 / d 3 x e - i(x - s) V 1/2 (S). (10) 

In order to verify this result and investigate the importance of the sign prob- 
lem, we have computed the discretised Wigner function (Q) by direct Monte Carlo 
sampling of the histogram for n. Each path comprises L independent points 
rii, i = 1 . . . L, taken from a uniform distribution on the sphere. Since the prod- 
uct of overlaps in Eq. ^ is complex for L > 2, a phase problem arises. For the 
results shown in Fig. [l , 10 7 independent paths were generated for each L from 2 
to 15. Further analysis shows the numerical distributions converging to the Wigner 
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Figure 1. Discretised Wigner function for spin 1/2 and 1 with L = 2 ... 15 time slices and 10 7 Monte 
Carlo steps. 



function (|3|) according to Eq. (|10(). We also see how the phase fluctuations lead to 
numerical instabilities. The correct distribution would be obtained if the the number 
of Monte Carlo steps were taken to infinity before the number of time slices. It is 
clear from the figure that 10 7 steps are insufficient for convergence for spin 1; results 
for higher spin show still slower convergence. 

We have thus demonstrated the emergence of a non-positive quasiprobability 
(the spin Wigner function) from the Berry phases of the integral. Numerical calcula- 
tions are however rapidly overwhelmed by the sign problem. 
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